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Abstract. In this article we study the fractal Navier-Stokes equations by using stochastic La- 
grangian particle path approach in Constantin and Iyer [6 |. More precisely, a stochastic rep- 
resentation for the fractal Navier-Stokes equations is given in terms of stochastic differential 
' equations driven by Levy processes. Basing on this representation, a self-contained proof for the 

existence of local unique solution for the fractal Navier-Stokes equation with initial data in W' '' 
is provided, and in the case of two dimensions or large viscosity, the existence of global solution 
is also obtained. In order to obtain the global existence in any dimensions for large viscosity, the 
gradient estimates for Levy processes with time dependent and discontinuous drifts is proved. 



o 

(N 



o 



1. Introduction 

Consider the following incompressible fractal or generalized Navier-Stokes equation in 
(abbreviated as FNSE): 



^ I r dfU = JLu - {u ■ V)u + Vp, f > 0, 

V • M = 0, m(0) = Mo, i^-^) 



where m = (m\ • • • , u'^f denotes the column vector of velocity field, p is the pressure, X is the 
generator of a Levy process given by 

£u(x) = I (u(x + y)- u{x) - l|v|<iCv • V)M(x))v(dj), (1.2) 

where v is a Levy measure on R'', i.e., it satisfies that v{0} = and 

r 1 A MVdy) < +CX,. 
- - - Jr'' 

When v(dj) = dy/|y|''^" with a e (0,2), X = -Cq,(-A)"^^^ is the usual fractional Laplacian 
operator by multiplying a constant. 

As a simplified model of equation (11.11) . the following fractal Burgers equation has been 
studied by Biler, Funaki and Woyczynski [3] and Kiselev, Nazarov and Schterenberg [[TSl . 

dtu = -{-^f''^u - (u ■ V)u, t > 0, m(0) = Mo- 

As for generalized Navier-Stokes equation (|1.11 ). when X = -(-A)'^''^, it has been studied by 
Wu ll25l in Besov spaces by using purely analytic argument. The main feature of such fractal 
equations is that operator X given by (11.21) is non-local. Recently, there are increasing interests 
for studying such fractal equations or fractional dissipative equations since they naturally appear 
in hydrodynamics, statistcal mechanics, physiology, certain combustion models, and so on (cf. 
ED 13 Ea, etc.). 
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The aim of this paper is to study equation (11.11 ) by using a stochastic Lagrangian particle 
trajectories approach following [|6] |27]. More precisely, Constantin and Iyer [[6l gave the fol- 
lowing elegant stochastic representation for the regularity solution u of Navier-Stokes equation 
(corresponding to X = vA in (|l.ll) ): 

( C 

Xrix) = x+ I u,(X,{x))ds + y/lvBt, t > 0, 

Jo (1.3) 

u, = FE[iVX')(uooX;')], 

where P denotes the Leray-Hodge projection onto divergence free vector fields, B, is a Brownian 
motion, and X~\x) denotes the inverse of x Xt{x). Basing on this representation, a self- 
contained proof of the existence of local smooth solutions in Holder space was given by Iyer 
|fT3l . Later on, by reversing the time variable, in a previous work [|271 , we considered the 
following stochastic representation: 

Xi ,(x) = x+ f UriXt r{x))dr + V2y(5, - B,), t^s<0, 

Jt ' (1.4) 

M, = PE[(V%,o)(woo^r,o)], 

and a self-contained proof of the existence of local smooth solutions in Sobolev space is also 
obtained. Moreover, the global solution for large viscosity is proven by using Bismut formula. 

Naturally, if one replaces the Brownian motion in (11.41 ) by a general Levy process Lf, then it 
is expected that the corresponding solution u will solve the following backward fractal Navier- 
Stokes equation: 

dtu + £u + iu-V)u + Vp = 0, t< 0, 

(1 .5) 



V • M = 0, m(0) = Uq, 

where X is the Levy generator of the Levy process L,. In general, it seems hard to solve the 
above fractal Navier-Stokes equation by using purely analytic tools. However, stochastic system 
(11.41 ) is easier to be dealt with if one replaces 5, by a general process and only considers the 
local smooth solutions. In fact, it is easy to obtain the existence of local smooth solutions for 
stochastic system (11.41) , and a global smooth solution in two dimensional case by the same 
arguments as in [T21I2S1. This will be given in Section 2. 

On the other hand, if we only assume that the initial data belongs to the first order Sobolev 
space W'''', it seems not so easy to construct a local solution in W'''' for stochastic system (11.41 ). 
A clear difficulty is to obtain the differentiability of the solution flow x i-^ ^f,o(.^)- Although 
one can solve the following equation with divergence free vector field u e L;'^^^.((-oo, 0]; W^ '') 
by using DiPema-Lions' theory (cf. Iil0n 7il): 



Xt,s(x) = X + J u,-(Xt^r(x))dr + (Lv - L,), 

it is only known that x i-^ Xf o(x) is approximately differentiable (cf. [HT]). This difficulty will 
be overcomed by using Krylov's estimate for jump-diffusion processes and the regularizing 
effect of Levy process if is non-degenerated in some sense (see Condition (H)q, below). 
Thus, following Section 2, Section 3 will be devoted to the proof of the existence of local 
W^'^'-solutions. 

For proving the existence of global solutions in W' '' for large viscosity, we need some gradi- 
ent estimates for the SDE with Sobolev coefficients and driven by a Levy process. For this aim, 
we shall use some asymptotic estimates for the heat kernels of Levy processes due to Schilling, 
Sztonyk and Wang [20]. Our approach for the gradient estimates of SDEs is based on the a 
priori estimate for an integro-differential equation and the uniqueness of weak solutions. This 
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is the content of Section 4, and can be read independently. In Section 5, we prove the global 
well posedness for large viscosity. 

Lastly, we mention that other stochastic approaches for incompressible Navier-Stokes equa- 
tions can be found in the references [fT6l 51 El |8l |9i, etc.; and compared with the analytic argu- 
ments, one of the main advantages of representation (|1.41 ) is that it is convenient for numerical 
simulations (cf. [[T7l[T4ll '). This is in fact our main motivation for studying the stochastic repre- 
sentation of fractal Navier-Stokes equation (11.51 ). 

2. Stochastic representation for fractal Navier-Stokes equations 

We first fix some notations. Set No := {0} U N and ]R_ := (-oo,0]. For k e No, let = 
Cl(W'; R'') be the space of all /r-order continuously diff"erentiable vector fields on R'' with the 
norm 

k 

ll"llc* •= / . sup \V-'u(x)l < +0O, 

where denotes the j-order gradient, and | • | denotes the Euclidean norm. For ^ e No and 
j9 > 1, let W'^ '' = W'''P(W'; W') be the usual vector-valued Sobolev space on R'' with the norm 

k 

M\k.p ■= W^^'^Wp < +°°' 

where || • ||p is the usual L^'-norm in R''. 

Let us now recall some basic notions and facts about Levy processes on negative time axis. 
Let {L,),^Q be an R'^-valued Levy process on some probability space (Q, P), i.e., an R^'-valued 
stochastically continuous process with stationary independent increments and Lq = 0. By Levy- 
Khintchine's formula (cf. [2» p. 109, Corollary 2.4.20]), the characteristic function of L, is given 
by 

E(e'^-^') = exp I? ib-^ + ^'A^ + / ~ ^'^"^ + " ^luKilHd^) | =: (2.1) 

where iff(^) is a complex- valued function called the symbol of (L,),^q, and b 6 R'', A 6 R'' x R'' 
is a positive definite and symmetric matrix, v is a Levy measure on R''. Throughout this paper, 
we only consider the pure jump Levy process, and assume below that 

b = 0, A = 0. 

We remark that t ^ Lt admits a version still denoted by L, such that for almost all a», ? i-^ U{o^) 
is right continuous and has left limit, but, for fixed t, 

P{o) : Lt{oj) + L,^{oj)] = 0. 

Below, for ? < 5 < 0, define 

^t,, ■= o-{Lr - Lt : t < r ^ s). (2.2) 

Given u e C(R_; C^(R^; R'')), for x e W^, let Xt^,{x) be the unique solution of the following 
SDE: 

Xt,s(x) = ^ + Ur{Xt^rix))dr + (L, - Lt), t<s<Q. (2.3) 
It is easy to see that {X, .^(x), x e R'', t ^ s <,0} forms a stochastic C^-diffeomorfism flow, and 



VZ,,,(x) = I + VuriXM) ■ VX,,,(-«)dr, (2.4) 



where VX,, j(jc) = {djXl^{x))ij=i...^d, and dj denotes the partial derivative with respect to the 7-th 
component of x. 

Let N{t, Y) := 2ks<o lr(^i -Lg-), T G SiM!^) be the Poisson random point measure associated 
with (Lf)t^o. By Levy-Ito's decomposition (cf. [2, p. 108, Theorem 2.4.16]), one has 

Lf = I xN(t, dx) + I xN(t, dx), 

Juki J|a:|>i 

where N{t,dx) := N{t,dx) - tv(dx) is the compensated random martingale measure. For g 6 
CliW'; R), by Ito's formula (cf. [El p.226. Theorem 4.4.7]), we have 



g(Xr,,) = g(x) + [£g{X,r) + (Ur ■ V)g(X,,)]dr + Mf„ 



(2.5) 



where £ is the generator of (L,)f^o given by (11.21) , and (Mf^)^g[, q] is a square integrable martin- 
gale given by 

Ml : = 




[g{X,r-+y)-giX,r-m{dr,dy). 



We have 



Theorem 2.1. Let u 6 C(R_; C^(R'';]R'')) anJ Xt^,{x) be the solution of SDE (Q. For (p 6 



C2(R^; R'^) and c e C(R-; Cl{W^\ R)), 



/i((^) := E 



exp 



(x))drU(X,o(x)) 



and 



w,{x):=m%fl{x)-ip{X,fl{x))l 
Then h,w e C^(R_; C^(R'';R'')) respectively and uniquely solve the following partial integro 
differential equations (PIDE): 

dth, + Lh, + {ut ■ V)h, + c,ht = 0, ho{x) = (p(x), 

and 

dfWt + £wt + (ut ■ V)w, + (V^Ut)Wt = 0, wo(x) = (p(x). 
Proof Fix ? < 0. For ^ 6 C^(R'') and 5 > 0, by Ito's formula (see (|23])), we have 



(2.6) 
(2.7) 



exp 



JCr{Xt-S,r( 
t-S 



^sAx))dr\ g{X,^sM)) 



exp c,(X,_,,,(x))dr| [Lg{X,.sAx)) + (i<. ■ V)^(X,_,,,(x))]d^ 



+ . 



exp 



t-S 



Cr{Xt^S,r{x)W \ C,{Xt-sAx))g{Xt-sM)) 



By the stochastic continuity of f 1-^ L,, from equation (12.31) . it is easy to prove that {t, s) X,^six) 
is also stochastically continuous. Thus, since {s, x) ^ £.g{x) + {Us-V)g{x) + Cs{x)g{x) is bounded 
and continuous, by the dominated convergence theorem, we have 



exp 

(5^0 



Noticing that 



J^c,(X,.,,,(x))dr|g(X,.,,,(x)) 
£g{x) + {u, ■ V)g{x) + Ct{x)g{x). 



Xt-sfl{x) — X,Q o Xt-s^t(x), 
4 



-8(x) 



(2.8) 



by Markov property, we have 



hf-s = E 



exp 
exp 
exp 



CriX,-s,r)dr^ exp |^ c,(Xj,, o X,_5,,)dr| (p(X,,Q o Xtsj) 

Cr(X,^S,r)dr^E e^pj^ Cr(Xt,rOXt-Sj)dr^(p(XtfiOX,^sj) 

I CriXf^gf) 
Jt-S 



Thus, by (|2.8I) . we obtain 

\{ht{x) - ht-six)) -[£ht(x) + {ut ■ V)h,(x) + Ct(,x)htix)]. 



Since the limit is a continuous function of (t, x), it follows that for each x, t h,{x) is differen- 
tiable and equation (12.61 ) is obtained. 
As for (12.71 ), observing that 

VX,.,,o(x) = (VX,o) o X,_g,{x) ■ VX,.g,{x), 
by Markov property again, we have 

w,.s{x) = E[V%.,,o(^) • ipiXr.sdx))^ = E[V%_,,(x) • w,{X,.sAm. 
Hence, by (|Z8]) and (l24l) . we have 

-(wt{x) - wM) = -lE[w,{Xt.sAx)) - wr(x)] - ^E[(V%_,,, - I) • WtiX,.s,(x))] 

o 



6^0 



-[£wt + (u, ■ V)w,]ix) - (V'ut ■ Wt)ix). 



Equation (12.71) is thus obtained. 

We now prove the uniqueness. Here we adopt the duality argument. Let Xt^x) solve the 
following SDE: 



Xt^x) = X 



- Ur{Xt,r 



Fix e C^(W; W) and T <0. Forte [T, 0], define 



ht{x) :-- 



exp 



JT't 



)(Xr-,.W)drU(Zr-,W) 



By the above proof, it follows that h, e L^W') n C^CR^) solves the following PIDE: 

dfit = JL*ht- {ut ■ V)h, + {ct - di\Ut)ht 
= -Cht - div(M, (g) hf) + Cfht 
subject to hj{x) = (pix), where X* is the dual operator of X and given by 

rgix) = r {gix -y)- g{x) + (y ■ V)^Wl|,|,i]y(d3;). 

Now, let h e C^R-; C^(]R^;]R^)) solve (IZ6l) with ho(x) = 0. Then, by the integration by parts 
formula, we have 

dtiht, ht) = -{£ht + {ut ■ V)ht + c,ht, ht) + {ht, -CK - di\{u, ® h,) + cA) = 0, 
where {ht,ht) = j^j{ht{x),ht{x))f^ddx. Since {ho, ho) = 0, it is immediate that {hT^hj) = 



QiT, (p) = 0, which then gives hrix) = by the arbitrariness of 6 Cq(W; W). 
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Remark 2.2. If one assumes that u e L/^^,(R-; C^(R^; W')) and c e Ll^.(R.;CliW'; R)), then the 
conclusions of Theorem \2.1\ still hold if one replaces equations A2.6\) and (EZl) by 

htix) = ip{x) + [£hs{x) + {Us ■ V)h,{x) + c,{x)h,{x)]ds, 

and 

Wt(x) = tp(x) + I [JLws(x) + (Us • V)ws(x) + (V^Us)Ws(x)ds. 



Using Theorem 12. 1[ we have the following representation for the solution of fractal Navier- 
Stokes equation as in [27, Theorem 2.3]. 

Theorem 2.3. Let u e C(R_; C^(R^; R'')) be divergence free. Then, u is a solution of fractal 
Navier-Stokes equation if and only ifu solves the following stochastic system: 

Xt^,(x) = x+ I Ur(X,,.(,x))dr + (L, - L,), t < s 

ii, = PE[V%,o-(woo«], 

where L is a Levy process with generator X, and P is the Leray-Hodge projection onto diver- 
gence free vector fields. 



Along the completely same lines as in 11271 Theorems 3.8 and 4.2], one has the following 
result. The details are omitted. 

Theorem 2.4. For any k gMq and p > d, there exists a constant Co = Co(fc, p,d) > such that 
for any uo 6 W^^^'^CR'^; R"^) divergence free and T := -(Co||VMollyt+i,p)~\ there is a unique pair 
of(u, X) with u e C{[T, 0];W''^^'P) satisfying dZP]). Moreover, for any t e [T, 0], 

Wut\\k+i,p < Co||VMolk+i,p- 
In two dimensional case, one has that for all t e R_, 

\\u,\\k+2,p < Ci\\Uo\\k+2,p,k,p,t), 

where the constant C is increasing with respect to its first argument. In particular, there exists 
a unique global solution u 6 C(R_; W'^^^ '') to l\2.9\l in the two dimensional case. 

3. Existence of local solutions for FNSE with W^ '' initial data 

In the remaining sections, we mainly study equation (12.91 ) with uq 6 W^'''(R''; R''). For this 
aim, we assume that 

(H)q, Let be the Levy symbol given in (12.11) and satisfy that for some a e (0, 2), 

Re(<A(^)) X as 1^1 ^ oo, 

where a ^ b means that for some ci, C2 > 0, cib < a < C2b. 
Consider the following SDE: 

= ^ + Ur{XrAx))dr + v^'"(L, - L,), t < s ^ 0, (3.1) 

where m : R_ x R'' ^ R"^ is a bounded Borel measurable function, and with a little abuse of 
notations, v > denotes a positive constant which plays the viscosity role. 

We recall the following Krylov estimate for jump diffusion processes taken from [l 28] The- 
orem 3.7]. Although the theorem is given therein only for or-stable processes, it is clearly also 
valid for more general Levy processes considered in the present paper since the proof only 
depends on the gradient estimate (|4.5I ) below. 
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Theorem 3.1. Suppose that (H)q, holds with a e (1,2), and u is bounded by k. Let v(jc) solve 
equation di./l) . Ther 

independent of v 
U{{s,,S2]-LP{W')), 



equation di./l) . Then for any p > ^ and q > there exists a constant = Cuid, a, p, q, ij/) > 
independent of v > I and x e W' such that for a// -1 < f < 5i < ^2 ^ and f e 



f 



(3.2) 



friXrixMr 

where is increasing with respect to k, and ^t,s^ is defined by A2.2\l . 

The following lemma is taken from [[TTl p. 1, Lemma 1.1]. 

Lemma 3.2. Fix t < 0. Let {jS('^)}.ve[f,o] be a nonnegative measurable {^t,s)-adapted process. 
Assume that for a// ? < < 52 < 0, 



ns2 


\ 




\ /3{r)dr 




<p{Si,S2), 









where p(si, S2) is a nonrandom interval function satisfying the following conditions: 

(i) p{si, S2) < p(.S3, S4) if(si, S2) c (S3, S4); 

(ii) \imsiQSup,^^^^^^^o,\s,-s2\<sP(^i^^2) = 0. 
Then for any 7 > 0, 

Eexply r /3ir)dr\ < 2^, 



where N eN is chosen being such that for any k = 0, - ■ ■ ,N - 1, 

p(-(k+l)\t\/N,-k\t\/N)<^. 

Let / be a locally integrable function on R.''. The Hardy-Littlewood maximal function is 
defined by 

A1/(jc) := sup — I f{x + y)diy, 

r>0 \tir\ JBr 

where B,. := {3; 6 R.'' : |3^| < r} and 15^1 is the volume of Br. 

We recall the following well known results (cf. [[TSl Appdenix] and [|22l p. 5, Theorem 1]). 

Lemma 3.3. (i) For any f e W^''', there exist Cd > and a null set E such that for all x,y iE, 

\fix) - f{y)\ < Q(M|V/|(x) + M^f\{y))\x - y\. (3.3) 

{ii) For any p > I, there exists a constant Cd,p > such that for any f e LP(W^), 

IIM/llp < CdJfWp. (3.4) 

Using the above three tools, we can derive the following important estimates for later use. 

Lemma 3.4. Suppose that (H)q, holds with a e (1,2), and p > ^. For any U > 0, there 
exists a time T = T{U) 6 [-1,0) independent of v > \ such that for any divergence free 
u e L'^([T, 0];W^'P{W'; R^)) with 

sup < U, (3.5) 

te[T,Q] 

the unique solution Xf ^(a:) to SDE ({3.1\l belongs to H^^i W]^^ with respect to x, and preserves the 
volume, and satisfies that for any y ^ 1 and some C = C(T, y, U) > 0, 

sup supE|VXf,o(x)r < C, (3.6) 

t&Vrm .veR'' 



and 



sup supE|VX,,oWr <2. (3.7) 

«e[r,0] xeWi 



Moreover, for any ip e W^'''(W^; W'), if we define 

w, := PE(V%,o ■ ° ^r,o)), (3.8) 

then we C([r,0]; W^-^) and 

dtw, = PE[V%,o • (V^ - VV) o X,o ■ diXrfll (3.9) 

Proof Under (l33]), it has been proven in fW. Theorem 1.1] (see also (TT\) that SDE (IXTl) 
admits a unique strong solution X, i(x) for each x 6 W'. Since w is divergence free, x i-^ Xf,A('<^) 
preserves the volume. Let uf(x) := Ut * Pe{x) be the mollifying approximation of u, where 
(Pe)£e(o,i) is a family of moUifiers. Let Xl^{x) solve the following SDE 

XlXx) = + J M^(X,';,(x))dr + vi/"(L, - L,), t<s<0. 

Then 

vz;;Xx) = I + J Vm^(z;;,(x)) • vz^,(x)dr, 

and 

\vxi{x)\ < 1 + J |Vii^(x,';,(x))| • |vx,';,(.^)|dr, 

where | • | denotes the Hilbert-Schmidt norm for a matrix. By Gronwall's inequality, 

\VXl(x)\ < exp U |V<(X;;,(x))|dr| . (3.10) 
By Theorem 13. 1[ one has that for any q > and all -1 < f < 5i S2 <:0, 



^ ^II"IIl~([,.0];W1-P)'I""^''([-*I.'*2];W'-P) 

where the second inequality is due to the Sobolev embedding relation W' '' c L°° provided that 
p > d. Hence, by Lemma [X2l for any y > I, 

sup supEexp jy r |VM^(Xf^(;c))|dr| < +c», (3.11) 

£6(0,1) .rf [ Jt ' ) 

and one can choose T e [-1,0) depending on U such that for all t 6 [T, 0), 



sup sup E exp ^ 4 f \Vu^^(XlXx))\dr \^2. (3.12) 



Let v(;c) solve the following linear random ODE: 



(Claim): For any 6 e [1,2), 



lim sup supE|Xf,(;c)-X,,,(x)|^ = 0, (3.13) 

t<s-t,se[T,0] xeBJi 



and for each t e [T, 0) and x e W', 



lim E| VXfo(x) - A,,o(x)r = 0. (3.14) 



We first prove (I3.131 ). By (13.31 ). we have 

< Q J iM\V4\(Xl(x)) + MWu^XrAxMXlix) - XrAx)\dr 

+ J lu'.iX.Ax)) - UriXUxWr, 
which yields by Gronwall's inequality that 

\Xl(x) - X,Ax)\ < exp |q J (M\Vu%Xlix)) + M|VM^|(Z,,(x)))dr 

X J KiXUx)) - UriXrAxmr. 
As in estimating (13.111) . one has that for any y > I, 



sup sup E exp <j y ( M\Vu^\(XlXx))dr \ < +oo 

£e(0,l)jc6R'' 



and 



sup supEexp<{7 ( M\Vu%XtAx))dr\ < +oo. 

ee(0,l).veR'' 



Hence, by Holder's inequality and Theorem [XT] again, we have for any 6 e [1,2) and q > — ^ 
E\Xl{x)-X,Ax)f <{Eexpl^^J^ M\Vu%XtAx))dr^^ 

X |Eexp 1^ J M|Vi<^|(X,,W))dr 



pa— 2d ' 

(2-a)/4 



X 



;J K - u/{x,Ax))d, 



SI2 



,£ .,|2||<5/2 



= Cll"'^ - ^^lli.,,,0];L.)' (3.15) 

where C is independent of e, x, t, s. Limit (|3.13|) now follows from the property of convolutions. 
As for (13.141) . we have 

\VXlix)-AUx)\ < J \M(.Kr(^)) - ^Ur(X,Ax))\ ' \VXl(x)\dr 
+ J |Vm,(X,,(-^))I • \VXfAx)-AAx)\dr, 

which then gives that 



\VXl,ix) - Ar,o(x)\ < J \Vu%Xl(x)) - Vur{X,Ax))\ ■ |VX,^,(.^)|dr 
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xexpjj |VM,-(X,,(x))|dr| . 
As above, using (13.101 ). (13.1 II ) and Holder's inequality, we have 

(^0 \'5/2 
E J \Vu';{Xl{x))-Vu,{X,,(x)fdr\ . (3.16) 

For fixed s' e (0, 1), by (I3.13D , we have 

Ej |V<(X,';,.(x)) - V<(X,,,-(x))pdr ^ as e ^ 0. (3.17) 

pa 



pa— 2d ' 



By (|3.2I) . we have for (^r > 

sup E |V(< - M,)(X,';,(x))|2dr < C||V(m^' - m)II^2,([,o];l.) ^ as ^ 0. (3.18) 

ee(0,l) Jt 

Limit (13141) then follows by (13161) . (IXT71) and (l3J8l) . 

Using the above claims, by (13.101 ). (13.111) and (13.121) . one finds that X,^,{-) e Pi-^^iW'^^, and 
(13.61) and (13.71) hold. Moreover, if we define 

then, since P is a bounded linear operator in and x ^fo^-^) preserves the volume, it is easy 
to see that w"" e C{[T, 0] ; n^gNW^'''). Noting that 

= FVE(Xl, ■ (if o XI,)) = VEiVX.o ■ (V ° ^'o)) + PE(V'(^ o X^,) ■ Xf,), (3.19) 
we have 

= PE[VXl, ■ (V^ - VV) o XI, ■ diXl,]. 
Using limits (13.131) and (13.141) . formula (13.91) then follows, and meanwhile, 

lim sup li< - = 0. 

So, w e C([T, 0] ; W^-^). The proof is complete. □ 

The following lemma gives the continuous dependence of solutions to SDE (|3.1i) with respect 
to u. 

Lemma 3.5. Suppose that (H)^ holds with a e (1,2), and p > ^. For U > and T e [-1,0), 
let u,ue L-^dT, 0];W^ 'P(W'; W')) be divergence free with 

sup ||Mf||i,p < U, sup < U. 

Let X, X be the solutions of SDE ( 13. il) corresponding to u, u. Then for any 6 e \\,2), q > -^^^ 
and t 6 [r, 0], 

supE|X,o(x) -X,o(.^)r < CA\u - M|li.,([,o] L.)' (3.20) 

where Ci only depends on U, T, a, 6, p, q, d, ij/. Moreover, for any ip,(p & W^'''(W'; W^), let Wt 
and Wt be 
te{T,0), 



and Wt be defined as in ( 13.^1) corresponding to (if, X) and {(p, X), then for any q > ^^"^^ and 



\wt-Wt\\l'<C2\\cp-(p\\l' + C, J \\ur-ur\\l^dr, (3.21) 



where C2 ( resp. C3 ) only depends on U, T,a,p,q,d,i// ( resp. U,T,a,p,q,d,i//,\\ip\\ 1 pj. 
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Proof. Estimate (13.201 ) follows from the same calculations as in estimating (13.151) . Let us look 
at (|3.21l) . Using mollifying approximation and (13.191) . we have 

wt-wt = PE[V(X,o - ■ (V o X,o)] 

+ PE[V'X,o • {if o X,o - (p o X,o)] 
= PE[V'(^oX,o)-(X,o-Ao)] 
+ PE[V'X,o • (V o X,o - if o X,o)] 
+ PE[V%,-(^ot,o-^oX,o)]. 
By the boundedness of P in L'\ we have 

\\wr - mil,, < C||E[V(^ o X,,o) • (X,o - X,o)]||;, 

+ C\\E[V%^o-(voX,o-<potMp 
+ C\\E[V%fl ■ (ifi o X,^Q - ip o XrMp- 
By Holder's inequality, for any 6 e (p/(p - 1), 2) and some /3 = /3{6,p) > 1, we have 

< C||E[|VX,,o| • (MIV^KX.o) + M\V<p\{X,fl)) ■ K,o - X,M\p 

< C||||VX,ollz/^(n) • l|M|V^|(X,o) + M\Vip\{l,o)\\mn) ■ \\X,,o - Koh^aMp 

< C sup \\VXr,o{x)\\if^n) ■ sup \\X,,o(x) - Xt,o(x)\yiii) ■ l|M|V^|||p 

< C||m - uWiMm-LP) ■ IIV^IIp, 
where we have used that x i-^ X,^o(x), X,^o(x) preserve the volume. Similarly, we have 
l|E[V^(^ o X,o) • (X,o - XrMp < C\\u - u\\ W<p\\p, 

and 

||E[V't,0 ■ if o Xr,o - if o XrMp < CWif - ^\\p. 

The proof is thus complete. □ 

We are now in a position to prove the following main result of this section. 

Theorem 3.6. For any divergence free uq e W^'''(R'^; R'') with P > there exist a time T = 
TiWuoWi^p) < independent of V > 1 and a unique pairof(u,X) with u e C([r, 0], W^'^) solving 
the following stochastic system: 



Xt^,{x) = + ^ Ur{X,^rixMr + v^'"{L, - Li), t<s<0. 



(3.22) 

= PE[V%,o-("oo«]. 
Moreover, for some Co > 1 only depending on p, 

sup WUtWip < 3Co||Molll,p, 
te[T,0] 

and u satisfies di.5l) in a generalized sense, i.e., for all divergence free vector field (p e C^(R'^;R''), 

(ut, (f)) = <Mo, + J + ■ V)m,, 0)]d5, (3.23) 
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where 



£ 



igix) = f \g(x - v'l"y) - g{x) + {v'l"y ■ V)g(x)l\y\^Mdy). 



Proof. Set M°(jc) = Uo{x). For n 6 No, by Lemma [341 we recursively define 

Xlix) = x + J u%XlXxMr + v'l"{L,-L,), 

Let us estimate the W^'''-norm of u"^^ . First of all, by Holder's inequality, we have 



(3.24) 



< c|^^ \\VXl,{xX,^^^ ■ \\uo o Xl,{xX,^^^Ax 

< C sup \\VXl^{x)\y^a, ( r Eliio o Xl^ix^Ax 
= Csup||VZ«o(x)|L2(n)||iiollp, 



and by (IXgl) . 



llVi/r^l,, < ciiE[|vx;;oP • ivmq ox;;oi]||p < c sup iivx;;oW||^4(^)||VmoIIp, 

Hence, for some Co > 1 only depending on p. 



\\uT%,p < Coll + supiivx;;o(x)ii24(n))ll"olli,;, 



Now, taking U = 3Co||moIIi,p in Lemma 1X41 by induction method, there exists a time T = 
T(U) < independent of v > 1 such that for all n e No, 

sup ||ar^lli,P < U, sup supE|VX;'„(jc)|^ < 2, (3.25) 

(e[7-,0] te[Tfi] xeRd 

and for any y > I, 

sup sup supE|VX"q(x)P < +00. 

neNo te[Tfl] xeR^ 

Thus, by Lemma [331 one has that for all t 6 [T, 0], 



where the constant C is independent of n, m and t. By Gronwall's inequality, we get 



lim sup Wu'; - u'l'Wl^ = 0. 



Thus, by (IXISl) . there exists aue L~([r, 0]; W''') such that 



and 



lim sup \\u" - UtWp = 

te[Tfi] 



sup ||m,||i,p < C/. 

te[T,0] 



Let Z be the solution of SDE (13.11 ) corresponding to the above u. By taking limits for both sides 
of (13.241) and using Lemma [331 one finds that iu,X) solves (13.221 ). Moreover, the regularity 
of a 6 C([r, 0]; W^ '') follows from Lemma \TM and the uniqueness follows from Lemma [331 
As for equation (13.231 ). let Uq(x) := Uq * Pe{x) be the mollifying approximation of u^. Let 
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Mf(x) e C([r, 0]; nyt^oW'-P) be the corresponding solution of (I3.221 ). By Theorem 12.31 one has 
that for all divergence free vector field cp e C^iW^; R'^) 

«, cf>) = {ut„ 0) + J [«, X» + <« • V)<, 0)]d^. 
Taking limits e ^ and by Lemma |331 one then obtains equation (|3.23l) . □ 

4. Gradient estimates for Levy processes with drifts 

In this section we prove gradient estimates for SDE (|3.1I ). It will be used to obtain the global 
well posedness for stochastic system (13.221) with large viscosity v in the periodic case. This 
section can be read independently and has some interests in itself. 

Let 'P{W') be the space of all Borel probability measures on R"^. For e P(W'), define 

T,f(x) = f f(x + yMdy), f e Q(R^). 

It is clear that for any v 6 P{W^), 

T,TJ = T,.J, f 6 C,(r'), 
where // * v denotes the convolution between /i and v. 



For a > 0, let us consider the truncated symbol of in (iXT 



^a{0= r (l-e'^" + /^-xl|,,Ki)y(dx). 

Jo<|,v|s;a 

For any f < 0, by Bochner's theorem, there exists a unique probability measure e f (R'') 
such that 



r 



JVL'> 



e'^Xr(d^) = e'^-^^^ (4.1) 

and respectively, 6 ^(R'') such that 

e'f-%/dx) = (4.2) 

We recall the following result from [|20l Proposition 2.3]. 
Proposition 4.1. Assume that (H)q, holds with a e (0, 2). Then for any a,t > 0, 

and for any n 6 No, there exists a time to = to(n, d, a,if/) <0 and a constant C = C{to, n, d, a, ijj) > 
such that for any t e [to, 0), 

|V>|,|i/.,,WI < C\t\-^''^"^l\l + \tV'l"\x\r''~' ■ (4.3) 
For /3 > 0, let W^-'' := (/ - A)~^'^(LP) be the Bessel potential space with the norm: 

\\f\\f,,p:=\\(I-Af^f\\p. 

If /3 = k eN, W^'P is the same as the Sobolev space W''^'. Using Proposition l4.1[ we now derive 
the following useful result. 

Theorem 4.2. Assume that (H)^, holds with a e (0, 2). For v > 0, define 

T;f(x):=E(f(v-^L, + x)). (4.4) 
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Then for any p 6 [1, oo] and m,n eN, there exists a constant C = C(p, m, n, d, if/) such that for 
all f 6 W'"-'', 

\Wf\Un,p < C[v{\t\ A Vv >0,\/t< 0. (4.5) 

Iff 6 (1, oo), the above estimate also holds for any nonnegative real numbers m,n. 

Proof. For v,a > 0, let fJ-l^pp-^^t ^ be defined as in (14.11) and (14.21) corresponding to the 

symbols ipaiv^^"^) and il/{v^'"^) - ij/aiv^^"^). It is easy to see that 

Ultidx) = v~" pa,,iv~" x)dx. (4.6) 

Set 

then for all ? e 

Vf = rifrij. (4.7) 

By (14.61) and Proposition 14. 1[ there exists a time t^ = ^0(^5 a, 1/^) < such that for all t e [?o, 0) 
and / e LP(R''), 

v"?:;;,/^ /(.^) = y-'S' f /(j)V'>|,|,/.,(v-^(j - .^))dy 

= v " f{x + y)S p\f^ua lv ''y)dy. 
Hence, for any p e [1, 00], by Minkowskii's inequality and (14. 3L we have 

liv%;,„/||;, < ''''^WfWpj^^ \^"p\,\"^^Ay'h)\dy 

< C(vk|)-(^"")/ll/||p r (l+iy\t\r'^''\y\r'-'dy 

Jr.'' 

= C{v\t\r'"\\f\\p r (1 + [yir'-'dy = C{v\t\r'"\\f\\p. (4.8) 

Jr'' 

Thus, by (|4.7I) and the -contraction of Tfl,,, we have 

\\^"^'"Vf\\p < C(v\t{r"'"\\t;^,^y"f\\p < C(y|?|)-"/l|V'"/||p, 

which yields that 

For general ? < it follows by the semigroup property of 7~/. If p e (l,oo) and m,n are 
nonnegative real numbers, it follows by interpolation theorem (cf. [|23l ). □ 

Let us consider the following PIDE: 

dth + £yh + iu-V)h = 0, r < 0, (4.9) 

subject to the final value 

hoix) = (fix), 

where 

£Mx) ■■= r + ''""y) - - (^'N ■ V)M(x)l|,|^ilv(d3;). 
By Duhamel's formula, one can write equation (14.91 ) as the following mild form: 

h,{x) = rjifix) + J rixius ■ v)/i.)(x)d5, (4. 10) 

where TJ is defined by (BT?! ). 
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We need the following simple lemma (cf. [|26l Lemma 5.1]). 

Lemma 4.3. Let z{t) be a nonnegative function defined on [-1,0), and satisfy that for some 
KuK2 > and/3,ye (0,1), 



z(t)<Ki\tr + K2 I 7-^ds, Vfe[-1,0). 
Then for any t e [-1,0) 

z(t) < CK,,/,,yKr\tr. 
We have the following useful estimate. 
Lemma 4.4. Assume that (H)q, holds with or 6 (1, 2) and 

u 6 L'"([-l, 0]; L'\R'^) + L'^iR'')) 

provided that p > Then for any ye + 1, c) and if e L^(R''), there exists a unique solution 
h € L'(r-1. 01: W^ P) to equation ii4.10\) such that for all t e [-1,0) andv> 1, 

\My,p<Ciiv\t\)'^^''M,„ (4.11) 

where Ci only depends on ■y,p,d,a,if/ and the norm o/ ||m||^,oo([_i o];LP(R'')+L~(Rrf))- In the case 
of u £ L°°([-1,0] X R.''), for any p e [l,oo] and (f e L''(W'), there exists a unique solution 
h 6 L^([-l, 0]; W^'P) to equation ([4.10\) such that for all t e [-1,0) and v > I, 

Wh,\\,<C2iv\t\)-'^''M,, (4.12) 

where C2 only depends on p, d, a, ij/ and the bound of u, and is increasing with respect to the 
bound ofu. 

Proof. We only prove the a priori estimates (14. lit and (14.121 ). As for the existence, it follows 
from a standard Picardi iteration argument. Assume that u = U1+U2 with u\ e L°°([-l, 0]; L^iR'')) 
and U2 6 L°°([-1,0];L°°(R'^)). Then, by (gj) and the Sobolev embedding W^-'' c provided 
(y - l)p > d, we have 

INIr,P < l|T,>llr,P + C J ivis-t))-h\(u,-V)h,\\pds 

< ||T;VII^,p + C I (s- 0^"(||mi,JIoo • l|V/i,||p + ||M2,.llp • l|V/i,|U)d5 

< c(v\t\)-^'"Mp + c J (s- 0""((ll"i,.IU + \\u2J\p) ■ \m,\Qds. 

By Lemma|43l we obtain (14.111) . 

In the case ofw 6 Li([-,0];L"(]R^)), one has 

Jo 
{v{s-t))--^\\(u,-V)h,\\pds 

<C(v\t\)-"''Mp + c J {s-t)-^\\u,\L-\\VK\\pds, 

which gives (14.121) by Lemma l43] again. □ 

Before proving the gradient estimates for SDE (13. It . we recall the notions of weak existence 
and uniqueness for SDE (13.11 ). Fix t < 0. Let D, be the space of all cadlag functions from 
[t,0] R''. We endow D, with the Skorohod metric so that D, is a Polish space. A weak 
solution of SDE (13.11) means that there exists a probability space {Q.,^,P) and two cadlag 
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processes (X,,.O.ve[/,o] and {Lt^s)se[t,o] defined on it such that s Lf^^ is a Levy process with 
symbol and 

X,^six) = x + UriXt^r(x))dr + v^'"L,^„ s e [t, 0], P - a.s. 

Such a solution will be denoted by {0.,^,P;{Xt^s)se[tm^i^t,s)se[t,o\)- Weak uniqueness means 
that for two weak solutions T^^, P^'">; (Xl[l)se[m, (Lfl)se[m), i = 1,2, the laws of ^ ^ 
in D, are the same for ?' = 1,2. 

Assume that for each x e R'', weak uniqueness holds for SDE (|3.1I) . Then for any bounded 
measurable function cp, the mapping 

X ^ E(p(Xtfi(x)) 

is well defined. Now, we can prove the following gradient estimate for x Eip{Xtfl{x)). 

Theorem 4.5. Assume that (H)q. holds with a e (1,2) and 

M eL"([-l,0];L^(R^) + L"(]R^)) (4.13) 

provided that p > We also assume that for each x e R^, weak uniqueness holds for SDE 
dJTD . Then for any y e + l,a), v > 1 and (p e LP(R^), 

l|E^(X,o(-))lly,;, < C,(v\t\r'"M„; (4.14) 
in the case ofu e L°°{[-1,0] x W'), for any p e[l,oo],v> 1 and (f e LP(W^), 

||VE^(Z,o(-))ll;, < C2(v\t\r''''Mp, (4.15) 
where Ci and C2 are the same as in Lemma \4~4\ 

Proof. Let uf(x) = u,*Pi;{x) be the mollifying approximation of u,, and Xf ^. be the corresponding 
solution of SDE dnj. For ^ e C~(R^), define 

h%x) := E^{Xl^{x)). 

By Theorem 12. 1[ u^ solves the following PIDE: 

dth' + £X + (a'^ • = 0, 

which is equivalent by Duhamel's formula that. 



h%x) = rjipix) + J TU(u: • v)/i:)(x)d^. 

Under (|4.13l) . it is well known that the laws of {Xf.(-'c)}ee(o,i) in D, is tight and any accumula- 
tion point is a weak solution of SDE (13.11) (cf. [28, Theorem 4. 1]). By the weak uniqueness, the 
whole sequence of XfXx) weakly converges to the weak solution Xt^.{x) in D,. In particular, for 

any x eW', 

lim /zf (jc) = lim E(piX^Jx)) = E(piX,^oix)). (4. 16) 

£—►0 £^0 ' 

By KWi . we have for all tp e C^(W^), 

\\E^iXl,i-)%^, = \\h%,^ < Civ\t\r^"M,. 

By (14.161) . we get (14.151) for any (f e C^(W^). For general (f e L'\W'), it follows by a standard 
approximation. 

In the case of m 6 L'"([-l, 0] x R''), by (|4.12D . for p = 00, we have 

\Ecp(Xl,(x)) - E<p(Xl„(y))\ = \h',{x) - h';(y)\ < C(y\t\)-""\\<p\Ux - y\. 
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Letting e — > yields that 

|E<^(X,oW) -E^(X,o(j))| < C{y\t\)-'l"MUx-y\, 

which then gives (14.151 ) for p = oo and ip e C^(R''). For general e U°(W^), it follows by a 
standard approximation. For p e [1, it is similar. □ 

5 . Existence of global solutions for FNSE with large viscosity 

In this section, we prove the global existence of W^'^'-solution for stochastic system (13.221 ) 
in the case of large viscosity and periodic boundary. Let T'' be the J-dimensional torus. Below, 
we shall work in the Sobolev space W''''(T^; R'') with vanishing mean denoted by Wq''(T''; R''). 
Thus, by Poincare's inequality, for any p > I, 

\\u\\p < C\\Vu\\p, 

and an equivalent norm in W''''(T''; R'') is thus given by 

INIi,;,-||Vm||;,. 

Below, we shall use ||Vm||p as the norm of Wj-'^CT^; R''). 

Theorem 5.1. For any divergence free uq e Wq'''(T''; R'') w?Y/z P > there exist v = v(||Molli,p) > 
1 sufficiently large and a unique pair of(u, X) with u e C(R_, W,)'^(T'^; R'')) solving stochastic 
system A3.22\i . 

Proof. By Theorem 13. 6 [ there exist a time T = T(\\Vuo\\p) e [-1,0) independent of v > 1 and a 
unique pair of (u, X) with u e C([T, 0], '') and 

\\Vu,\\p < 3CoI|VmoIIp 

solving stochastic system (13.221) . Our aim is to prove that for some = T^dlVwollp) 6 [7^,0) 
and V large enough, 

WVutMp < WVuoWp. (5.1) 

After proving this estimate, one can invoke the standard bootstrap argument to obtain the global 
solution. 

By (1X91) . we have 

diUt = PE[(V%,o - I) • Vuo o • diXt^o] + mE[uo o X,o] 

+ PE[V%,o • V'^o o • di(X,fl - X)] + PE[V(i<i, o X,o)]. 

Noticing that 

VXt^,(x) = J ^Ur(Xt,r(x)) ■ VZ,,,(x)dr, 
by Holder's inequality. Theorem ISTI and (13.61 ). we have for any ye (1' ff-* ^^'^ x eW', 

\VuriXtAx)tdr J \VXtAx)fdr 

\Vur{XtAx))\^^dr E J \VXt,r(xf^'dr 

where 7* = -pi and C is independent of x and t, and depends on IIVmqIIp- Hence, by Holder's 
inequality, 

||E[(V%,o-I)-VmooX,o-5A,o]|Ip 
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< C\t\\ 



< IlllVXf^o - MlHsi) ■ ll^"o ° Xtfl\\u'((i) ■ \\diXtfl\\i2p/(r-2)(^[i)\\p 

< C sup \\VXt,o{x) - I||L2(n) ■ sup ||o';Xf,o(x)||L2W(P-2)(fi) • llVaollp 

xeR'' xeR'' 

< C\t\ ■ WVuoWp. 

Similarly, 

||E[V%,o • V'uo o Xrfl ■ di(Xr,o - xMp < C\t\ ■ IIVMollp. 
On the other hand, by (14.151 ). we have 

\\diE[u,oX,^,]\\p + IIVEK oX,o]||;. < C(v\t\)-'/''\\u,\\p < C(v\t\r"''\\Vuo\\p. 

Combining the above calculations, we obtain that 

\\VUr\\p<C{\t\+(v\t\r'^"Wuo\\p, 

which then gives (15.11) by first letting t small enough and then v large enough. The proof is thus 
complete. □ 
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